Exact solution of one class of Maryland model 
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The Hamiltonian H of one-body Maryland model is defined as the sum of a linear unperturbed 
Hamiltonian Hq (Honn = nui) and the interaction V , which is a Toeplitz matrix. Maryland model 
with a doubly infinite Hilbert space are exactly solved. Special cases of one-body Maryland model 
include the original Maryland model (Phys. Rev. Lett. 49, 833 (1982) and Physica 10D, 369 (1984)), 
which describes a quantum kickied linear rotator and single band Bloch oscillations. Maryland 
model and single band Bloch oscillations are the same Hamiltonian in two different representations. 
A special case of many-body Maryland model is Luttinger model. 



PACS numbers: 05.45.Mt, 72.10.Bg 

I. EXACT SOLUTION 

The Hamiltonian of one-body Maryland model is de- 
fined as 



H = H a + V. 



(1) 



Both the unperturbed Hamiltonian Hq and the inter- 
action V are doubly infinite matrices; the indices n 
and m of ffo™ and V nm run from — oo to oo. Hq 
is a diagonal matrix; its diagonal matrix elements are 
{— oo, • • • , —2, —1, 0, 1, 2, • • • , oo} x uj(t). V is a Toeplitz 
matrix: V nm = V n - m . 

We want to solve Schrodinger equation 



(2) 



The unitary operator U (t, 0) dictates the dynamics of the 
system. 



4>{t) = u{t,o)i>(p). 

In the interaction picture, 



ih^U I {t,Q)=H I {t)U I {t,0). 



where 



Hi(t) = el/o'^Cl^e-i/o'^C')*', 
It can be verified, the matrix element of Hj(t) is 
Hr(t) —V pi("- m )/o w(t')dt' 

11 I\ L )nrn — v n—va^ 



(3) 
(4) 
(5) 
(0) 



Hj(t) is a Laurent matrix. See the Appendix A for a def- 
inition of Laurent matrix. From the theorem of Laurent 
matrix in the Appendix A, \Hj(t), Hj(t')] = 0. 
From the scattering theory, 



[//(£, 0) = Pexp 



Uo H{t ' )dt ' 



(7) 



where P is the time ordering operator. Since 
[H[(t), H[(t')\ = 0, we can remove the time ordering op- 
erator P here. 



Ui(t,0) = exp 



(8) 



The solution of the Maryland model Eq. (1) is expressed 
as the matrix exponential of a Laurent matrix. 

In the paper, we use Eq. (8) to solve two special cases 
of Maryland model, the original Maryland model p], [3, 
3] and single band Bloch oscillations [||. As the first 
Hamiltonian of the form in Eq. (1) was discussed in 
the model Hamiltonian in Eq. (1) is referred as 
Maryland model. Another term can be linear Toeplitz 
system, which emphasizes the matrix structure of the 
Hamiltonian and the role of the time ordering operator. 



II. CONNECTION WITH BLOCH 
OSCILLATIONS 

The physical explanation of the Hamiltonian in Eq. 
(1) can be the original Maryland model [3, 0, In 
the rotator representation [J 0, [H, H = —ifr-§g X ui(t), 
V = V(9,t) and \n) is -^^e m6 '■ In another physical 
explanation, |n) is seen as a site on an one dimensional 
lattice. V is treated as the kinetic energy which causes 
the electron hopping between different sites and Hq the 
potential energy of the electron in the linear electric field. 
The Hamiltonian of the quantum kicked rotator 



H 



1 8 2 
-ft 2 — 

2 ae 2 



k cos 9 S(t 

n=l 



(9) 



is just a lattice in a harmonic potential V(x) = ^h 2 x 2 . 
The kick strength k is the free diffusion time of the elec- 
tron. 

In the field of quantum chaos, the term "dynamic lo- 
calization" 0, Q means the absence of diffusion in the 
momentum space of the quantum kicked rotator when 
the kick frequency and the rotator frequency is incom- 
mensurable. Dunlap et al also used the term "dynamic 
localization" in the Bloch oscillations problem (See Sec- 
tion V) [4] . The result of the localization of the original 
Maryland model (See Section IV) is unexpected from the 
perspective of quantum chaos. But it is just Bloch oscil- 
lations. We use the term Bloch oscillations in a broader 
sense, while the ordinary meaning is an lattice electron 
in a time independent linear electric field. Even if the 



2 



hopping matrix elements between different sites are time 
dependent, under some condition, the electron still os- 
cillates on the lattice. The localization mechanisms of 
quantum kicked rotator [6| and the original Maryland 
model and single band Bloch oscillations are fundamen- 
tally different. 



III. EXACT SOLUTIONS IN THE ROTATOR 
REPRESENTATION 

In Section I, we solve the Maryland model in the site 
representation, now we solve it in the rotator represen- 
tation. The Hamiltonian of Eq. (1) in the rotator repre- 
sentation is 



In a M dimensional space, the solution of the Maryland 
model 



A I 



H = ^] u m (t)p m + V(x 1 ,x 2 , ■ ■ ■ ,x M ,t) (17) 



m— 1 



IS 



exp | - V[xi + ux{t")dt",x 2 + uj 2 (t")dt h 



, XM + j LJ M (t")dt",t']dt' }, 



(18) 



H=-iu(t)h-^+V(6,t), 

where the unperturbed Hamiltonian 

d 

H a = -iuj(t)h— = uj(t)p, 



(10) where x rn is the position and the momentum p m = 



(11) 



A special case of Eq. (17) and (18) is Luttinger model 
01 



JV 



M 



N M 



where the angular momentum operator;? = —ih-gg. Since H ~ 5^ Pi,n+ ^ P2,m + ^^^^V(xi, n X2, m )- (19) 



n— 1 m— 1 



n— 1 m— 1 



exp < — - / u(t ) dt'p 
ft Jo 



(12) 



is a translation operator, in the interaction picture, 
H T (t) = e iti aj( - t "> dt ' p V(6,t)e- : kti u( - t "> dt ' p 



r % N M i 

U (t, 0) = exp | - - (J2 Pl,n + ~P2,m)t | 
^ n— 1 m— 1 ^ 



V[ 9 + j u(t')dt',t 



(13) 



N M 



exp <^ - t(^2 V ( x hn ~ X 2 , m ))t \, 



Hi(t) is a function of position 9, so [Hj(t), Hi(t')] = 0. 
So the time ordering operator can also be removed here 
and 



0) =exp| - -J Vft + J u{t")dt",t']dt' 

The unitary operator in Schrodinger picture is 

U s (t,0) = e*ti Hodt Ui(t,0) 
= e —k / * <"(*') dt 'Px 

< 4 rt /■*' 



(14) 



(15) 



<>Xl>i - ~J V l + J L0(t")dt",t']dt' 



(20) 



where pi >n = ~^'o§[ — are position and momentum 
of the n-th of "1" particles (electrons) and x 2 ^ ni P2.m — 
—ih-Q^ — are position and momentum of the m-th of "2" 

particles (holes) and N and M are the total number of 
"1" and "2" particles respectively. The linear Tocplitz 
structure of Hamiltonian in Eq. (19) may be the origin of 
anomalous properties of Luttinger model and Luttinger 
liquid compared with Fermi liquid @, B Q ■ 



IV. THE ORIGINAL MARYLAND MODEL 

The Hamiltonian of the original Maryland model or 
quantum kicked linear rotator (QKLR) is 



Since 



d 



{H!(t)) nm = i- y * v {° + J q u{t')dt'^ m -^ e de 

— V f > i (" _m )(/o ""(*') dt ') 

(16) 

where V nm = V n . m = ± V{6, t)e^ m -^ e dO. The 
equivalence between the methods in Section I and this 
Section can be easily verified. 



(21) 



where 9 is the angle of the rotator, r the kick period 
and fc the kick strength, p = —ih-§g is the angular mo- 
mentum and V{9) — k cos 9 [Jl1,[1|- QKLR is unphysical 
because there is not a ground state. The classical version 
of Eq. (21) is not chaotic. The phase space is filled with 
invariant curves |[. Does QKLR delocalize if 4- is ra- 
tional? Berry proposed the energy of the rotator grows 
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quadratically From the perspective of the exactly 
solved eigenstates of the Floquet operator [2, H, [(| , the 
extended eigenstates generally mean derealization. But 
we will prove in the section, QKLR always localizes ex- 
cept one case when r is an integer multiple of 2ir. We set 
h = 1 and restrict t to the domain (0,27r). The Floquet 
operator 



_ P -iV(e) e -ipr 



F = e 



(22) 



The matrix representation of the QKLR Hamiltonian is 



H 



-2 |/ 



\ 



V 



/ -1 If 

if o if 
if i §/ 

y 2 



(23) 



where / 
picture, 



f(t) = Y^=i ~ nT )- I n the interaction 



If ^ = |, the period of QKLR is q kicks (2irp). In 
quantum kicked rotator, F = exp(— ikcosO) exp(— ip 2 ^). 
When t = 2ir, F is the same with QKLR with r = ir. So 
the period of quantum kicked rotator with r = 2ix is 2r 

(4tt) a. 

Now we discuss the seemingly conflict between the uni- 
tary operator F N and eigenstates of the Floquet operator 
F calculated in P, 0, [J- We calculate eigenvalues and 
eigenstates using a method similar to Q . The eigenvalue 
equation of the Floquet operator F is 



F<j>{6) 



Since e lpr is a translation operator, 

e -tVW e -fPT^ _ e-ivW^e - r ) = A0(6»). 

We assume ^ = £ . From Eq. (31), 
e -ine-r)^ 9 _ 2t ) 



(30) 



(31) 



-iV(fl-2r) 



3r) 



A0(0-r); 
A0(0-2r); 



|e lt 



e 



V 







2 l 





/ 



(24) 



Hj(t) is a Laurent matrix. [Hj(ti), Hj(t 2 )] = 0. From 
Eq. (8), to calculate E7>(i,0) = exp[-i #/(t') di'], we 
need to calculate J Q Hj(t') dt' . When t = iVr, the only 
non-zero matrix elements of Jq^ 1 " H(t') dt' are 



if(t') dt' 



fce lT (l 

2 r 



fc sm(^ T ) i(N+1) T 
2 sin(f) 



(26) 



7e 



where 7 = t ^n^f and 5 = + iW 2 , and 



JVr 



ff(t') dt' 



7e 



,-i<5 



(27) 



n,n+l 



Since 7 < 2 s in(V/2) ' ^(*>0) ^ s a (almost) band matrix. 
The rotator localizes, whether is rational or irrational 
except when t is an integer multiple of 2tt. But the 
localization mechanism is different from quantum kicked 
rotator [(| as we discussed in Section II. 
From the Appendix B, 



U(NT,0) Inm = e-^ m -^ 5 i m - n J, 
In Schrodinger picture, 

U(NT, 0) S nm 



,(27)- (28) 



< -- e l ( m -"^z m ~"J„_ m (27). (29) 



-inNr , 



e -iV{e-{ q -x)T) (j) (Q_ qT } = Xtf>(6-(q-l) T ). (32) 
From Eq. (31), (32) and <f>(d - qr) = 



q-l 



exp -nr)) = A". 



(33) 



n=0 



9-1 9-1 

^ 7(6» - nr) = ^2 k cos(6 - nr) 

n=0 n=0 

2 sm(^) 

= 0. 

Another derivation of Eq. (34) is 

q-l .-q-l 

k cos((9 — nr) = kRe exp(i(8 — nr)) 

q-l 



(34) 



0, 



exp(i6*) exp(z(— nr)) 



(35) 



where i?e is the real part of a complex number. So 

A" = 1. (36) 



A = 1, e 9 , e « , • • • , e 



(37) 



The Floquet operator F 1 has only g eigenvalues. F for ra- 
tional 27 has infinite degenerate point spectra. Grempel 
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et al thought quasienergy bands have a finite width 
Now we find eigenstates. Let's restrict 8 to the domain 
[0, 2s], Define <p(6) in the domain [0, as an arbitrary 
function and in other domains 



4>(e-r) = Xe tV( - e U(0); 
<f>{9 - 2r) = Xe lV ^ e - T U{0 - r); 

- (q - 1)t) = Ae iy ( e -(«- 2 ) r V(^ - (? - 2)r). 



(38) 



Eq. (38) gives the eigenstates of F. 

When is irrational, the eigenvalues are e llT , where 
I is an arbitrary integer, and the eigenstates are local- 
ized 0,0,13. From Eq. (38), we can construct localized 
and extended eigenstates. There are no real conflict be- 
tween the unitary operator in Eq. (28j^and (29) and the 
extended eigenstates calculated in 



(28) and 

tail. 



V. TIME DEPENDENT H 

The time dependent Hq considered by Dunlap et al 



H(t)=T (\m){m + l\ + \m + l){m\) 

m— — oo 

oo 

+ E(t) V" m\m)(m\, 



(39) 



m— — oo 



where T is the nearest-neighbor coupling and E(t) is the 
time dependent linear electric potential. Dunlap et al 
gave the analytic solution of the above Hamiltonian Q 
and found the time dependent field H generally destroys 
Bloch oscillations.. Here we treat Eq. (39) as a Maryland 
model or a linear Toeplitz system. 



(40) 



#(*) = E E T n (t)\m + n)(m\ 

n m=—oo 

oo 

+ E(t) ^ m\m){m\. 

m— — oo 

In the interaction picture, 

oo 

H *(t)=J2 E T n (t)e*ft B W*\m + n)(m\. (41) 



n m— — oo 



Hi(t) is a Laurent matrix. U(t,0)i = 

exp[— i J Q Hj(t') dt']. In the simple case of Q, 
Ti(t) = T_i(t) = T and E{t) = -Esmt. In N 
periods from to 2Nir, for n = 1, 

/ r^t'yti E(t")dt" di < 

JO 

~1 ATtt 

/*' -saint" dt" dt f = 2 7 rNTJ {E)t 



P-2N-K 

= T e 



-iE. 



(42) 



and for n = — 1, 



(43) 

In the derivation of Eq. (42) and (43), we used the for- 
mula 



iz cos 6 



J2 Jn{z)i n e m \ 



(44) 



n— — oo 



of [10|, From the Appendix B, the matrix elements 
of the unitary operator in the interaction picture is 

U(t, 0) Jnm = e^ n - n ^ E i m ~ n J n - m (^NTJ {E)). (45) 

Since 

e i^»H (t')df = lj (4g) 

the unitary operator in Schrodinger picture 

U{2Ntt 1 Q) s = U(2Ntt, 0)/. (47) 

From Eq. (46) and (47), the electron will diffuse away 
except when Jq(E) = 0. In the case of Jq(E) = 0, the 
electron will not delocalize. This is referred as dynamic 
localization by Dunlap et al 0]. 

VI. CONCLUSION AND DISCUSSION 

In summary, the Maryland model with doubly infinite 
Hilbert space is exactly solved. In the interaction repre- 
sentation, the unitary operator is the matrix exponential 
of a Laurent matrix. It is the special structure renders 
the Hamiltonian solvable. We think the solution can be 
generalized to a more general structure of Hamiltonian. 
We give the correct solution of the original Maryland 
model, concerning the resonant cases. Compared with 
Dunlap et aVs method [1| to solve the lattice electron in a 
time dependent linear electric field, our method based on 
the Maryland model is physically appealing and simpler. 
Further work should generalize the structure of linear 
Toeplitz system and remove the requirement of doubly 
infinite Hilbert space. 



APPENDIX A: 



A THEOREM OF LAURENT 
MATRIX 



"Doubly infinite dimensional Toeplitz matrix" is re- 
ferred as Laurent matrix in the mathematical literature. 
Two Laurent matrices commute. In the appendix, we 
give a proof (it may exist in another place) of the com- 
mutativity. A Laurent matrix A is defined as A nm — 
A n +i,m+i = A n ^ m , where i is an arbitrary integer and 
n, m run from — oo to oo. C = AB. C nm = A n iBi m . 
D = BA. D„ 



A 



5 



only M n+ i^ n and M n ^ n+ i are not zero. 



Cnm D nm ^ A n iB\ rn ^ ^ B n pAp m 

I p 

= ^ ] A n -iBi- m — ^ ' B n —pAp- m (Al) 
i P 

I p 

Given a I, there is a, p = n + m — I, which satisfies 
n — I = p — m and I — m = n — p, and vice versa. 
So J^i A n -iBi- m and ^ p i p _ m B n ^ p contain the same 
terms. C nm - D nm = 0. C = D. [A, B] = 0. Note 
the condition of doubly infinite dimension is necessary to 
ensure p = n + m — I is always the index of an existent 
matrix element. 

In fact, the above theorem is trivial because every Lau- 
rent matrix A is a function of the position 9 (a multipli- 
cation operator) in the position representation. 

oo 

A(0)= A n e- ine . (A2) 

n=— oo 

Two functions of position commute with each other. 



APPENDIX B: MATRIX EXPONENTIAL OF A 
LAURENT MATRIX 

We now calculate the matrix exponential e~ %M of a 
simple Laurent Hcrmitian matrix M. M is bidiagonal; 



M n+ i tfl = ~fe 



(B2) 



If we choose the basis of M as —l=e me , in the rotator 

V2tt 



M„, n+ i - 1 e- t \ 

where 7 and 8 are real numbers 

If we choose 
representation. 

M{6) = je l5 e i9 + ie - lS e- w = 2 7 cos(<9 + 8). (B3) 

The matrix element of e~ lM is 

(e~ iM ) 

1 



271 



= _L / e -iM(e) e i(m-n)9 de 



2tt 



2vr J Q 



2tt 



-i2-f cos(0-\-S) ^i(m — n)0 



-i2^f cos (0-\-S) ^i(m — n)(9-\-S) ^—i(m—ri)5 



(B4) 



= e -(™-»)V-"4- ra (-27) 

= e-^ m -^ s i n - m {-l) n - m J n ^ m (2j) 

= e- 4 (™-")^ m -"J n _ m (2 7 ), 
where J is the Bessel function of the first kind. In the 
derivation of Eq. (B4) , we used Eq. (44) and 



J n (z) = (-l) n J n (-z) 



(B5) 



of mini. 
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